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We show that in a ring-shaped Luttinger Liquid (LL) in 
contact with an electron reservoir — the chemical potential 
in the ring being controlled by a gate voltage Vg — both 
the average ring charge and the persistent current in the ring 
are step-like functions of Vg at low temperatures. The step 
positions are determined by the LL parameter a — vf/s, 
which therefore can be directly measured. We study electrons 
both with and without spin, taking into account long-range 
interactions in the ring. 



I. INTRODUCTION 



In particular, under conditions of strong transverse con- 
finement, when the charge motion in a nominally 2D 
quantum channel becomes essentially one-dimensional, 
one can expect to be able to realize the Luttinger liquid 
(LL) regime. Experimental studies of transport proper- 
ties of laterally confined 2DEGs do not provide unam- 
biguous proof of LL-like behavior of charges in quantum 
wires.Q Therefore it seems interesting and important to 
examine thermodynamic properties of quantum dots and 
wires and to suggest experiments which could reveal a 
non-Fermi-liquid character of the electron dynamics in 
systems of reduced dimensionality. 
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Quantum transport of charge in mesoscopic systems 
has been intensively studied — both theoretically and 
experimentally — during the last decade. Two kinds of 
oscillatory phenomena have frequently been in the focus 
of interest: 

i. Coulomb blockade of tunneling through small ca- 
pacitance metallic grains and the oscillatory behav- 
ior of the conductance upon change of gate voltage 
(for a review see Ref. ^ 

ii. Aharonov-Bohm oscillations of thermodynamic 
(persistent current) and transport (conductance) 
properties of mesoscopic rings in weak magnetic 
fields (see, e.g., the review articles Refs. |[||). 

Recently, both these effects were observed in semicon- 
ductor nanostructures.ffl Unlike in the case of metals, 
semiconductor heterostructures enable one to study the 
properties of electron systems of reduced (2D and ID) 
dimensionality. In nanostructures created by laterally 
constraining the motion of the 2D electron gas (2DEG) 
formed in an AlGaAs/GaAs heterostructure, both the 
characteristics of the system under measurement (elec- 
tron density, geometric form and size) and the the na- 
ture of the leads connected to it can be controlled. This 
is achieved by varying (i) the voltage, Vg, on a gate that 
couples capacitively to the net charge in the device, and 
(ii) by adjusting voltages on additional (split gate) elec- 
trodes that are needed for creating the lateral confine- 
ment. Hence, one is able to vary the parameters of the 
electronic device at hand and_its interaction with electron 
reservoirs over a wide range.Q 

In systems of reduced dimensionality electron-electron 
correlations play a significant role and the naive picture 
of noninteracting (or weakly interacting) electrons mov- 
ing in an external electrostatic potential does not apply. 
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FIG. 1. Sketch of the ring-shaped conductor discussed 
in the text. It encloses a magnetic flux <1> and is created in 
the 2D electron gas of a gated semiconductor heterostructure. 
The ring is in weak contact with an electron reservoir, which 
implies that the number of electrons in the ring is not fixed at 
flnite temperatures. The difference in electrostatic potential, 
Vg, between the ring-shaped conducting area and the rest of 
the 2D electron gas (the reservoir) is produced by a charged 
gate electrod. Because of the circular symmetry no electric 
field appears along the circumference of the ring. 

As a possible test of non- Fermi-liquid behavior of ID 
electrons we propose to study persistent current oscilla- 
tions and peculiarities of charge quantization in a quan- 
tum ring controlled by a gate voltage. The ring is as- 
sumed to be in contact with a reservoir of 2D electrons as 
illustrated in Fig. |l|. Therefore the number of ring parti- 
cles is not fixed at finite temperatures and interactions do 
affect transport properties even in an impurity-free ring. 
The voltage shift, Vg, between the conducting ring and 
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the rest of the 2D electron gas is produced by a charged 
gate electrod. In contrast to the situation dicussed in 
Ref. 1^ where the connection between ring and elec- 
tron reservoir is asymmetric, the exchange of electrons 
between ring and reservoir in our case takes place along 
the entire ring circumference. Hence, it is symmetric and 
does not lead to any electric field along the circumference 
of the ring. 

To calculate persistent currents and Coulomb blockade 
effects in a ring of strongly correlated ID electrons we use 
below the Luttinger liquid modeltJ both for spinless and 
spin- 1/2 fermions. 

We first show that the persistent current for electrons 
without spin oscillates as a function of gate voltage, Vg, 
the ring periodically exhibiting a dia- or paramagnetic 
response to the magnetic flux. The properties of these 
oscillations are sensitive to the value of the correlation 
parameter a = vp/s of the LL {vf is the Fermi velocity, 
s is the plasmon velocity). At low temperatures, when 
the amplitude of the persistent current attains its maxi- 
mum value Iq (/q = evp / L for an impurity- free ring, L is 
the ring circumference) the oscillations for a < 1 (repul- 
sive interaction) have the form of periodically repeated 
current bumps and dips. The width, i.e. the range in 
gate voltage of the bumps depends on the parameter a 
and the dimensionless flux $/$o (*i>o = hc/e is the flux 
quantum). In particular, at small flux values, (<&—!■ 0), 
the plateau width for paramagnetic response is maximal 
for strong repulsion (a ^ 1) and reduces to narrow spikes 
in the limit of noninteracting particles (a 1). For an 
attractive interaction (a > 1) the low-T oscillations of 
persistent current with gate voltage disappear. 

The above picture of current oscillations is correlated 
with the behavior of the average ring charge, Q, when 
the gate voltage is varied. For a repulsive interaction 
{a < 1) and low temperatures, the total charge in the 
ring is quantized in units of the elementary charge e. The 
function Q — Q{Vg) exhibits a typical behavior (Coulomb 
"staircase" ), which is a distinct signature of the Coulomb 
blockade. The step positions and their widths (for even 
or odd number of particles in the ring) depend on the cor- 
relation parameter a. For a stiff system (a ^ 1) all steps 
are of equal width. The width of the steps corresponding 
to an even number of spinless electrons is decreased when 
a is increased and collapses to a point when a — > 1. For 
an attractive interaction between spinless particles the 
mean ring charge is quantized in units of 2e and the os- 
cillations as a function of gate voltage vanish. 

At high temperatures, i.e. when T > T*^{a) — T*^ 
being the crossover temperature for current (j) or charge 
(c) oscillations — the oscillations of both persistent cur- 
rent and average charge in the ring become sinusoidal 
with an exponentially small amplitude, ~ exp(— T/T*). 
For current oscillations the crossover temperature is in- 
creased if the strength of the (repulsive) interaction is 
increased. In the limiting case of a stiff Wigner crystal 
(a 1) it has a maximum value T* = A/tt^ (for spin- 
less electrons; A = 2'KhvF/L is the energy level spacing 



at the Fermi level of non-interacting electrons at the rel- 
evant density) which coincides with, the crossover tem- 
perature for an isolated LL-ring£3ll3 Note that now and 
below vp always refer to the Fermi velocity of spinless 
electrons. The actual Fermi velocity for electrons with 
spin and the same density is therefore u_f/2. For the 
charge oscillations (as function of Vg) in the limit consid- 
ered (a — > 0) the crossover temperature goes to infinity 
as T* ~ l/a^. This means that at any finite tempera- 
ture a ring of strongly repulsive electrons exhibits sharp 
charge quantization. In the high-T region with strong 
attraction, > 3, we predict also that the period of 
the dominant contribution to the current oscillations is 
reduced by a factor of two. 

The above results correspond to the case of an in- 
teraction which is short-ranged. We show further 
that long-range Coulomb forces can be incorporated 
into our model if we replace the correlation parame- 
ter a of the LL by the effective "coupling" aeff{L) ~ 

a/ [1 + {2a^/TT)ie^/hvF)HL/X)]^^^, where A is the 
width of ID channel and L is the ring circumference 
(A<i). 

The magnetic field normally used in Aharonov-Bohm 
experiments is weak and the ring electrons are there- 
fore unpolarized even at low temperatures. To describe 
a real situation involving interacting fermions with spin 
one needs to understand how the electron spin affects 
the oscillation phenomena of interest here. We take the 
spin degrees of freedom into consideration by using the 
two-component LL model (spin- "up" and spin- "down" 
electrons of equal density). Our calculation scheme can 
be easily adopted to this model. For the purpose of 
comparing with the case of spinless electrons it is use- 
ful to distinguish between the cases with different par- 
ity of the number of spin- "up" (N^p) and spin- "down" 
i-^down) electrons (referring to the situation at zero tem- 
perature when the particle number is fixed). For a re- 
pulsive interaction there are three different cases: (i) 
(odd,odd) Nip = iVO„„„ = 2n + 1; (ii) (even,even) 
^up — ^down ~ 2n; (ui) (even, odd) or (odd, even) for 
N°p = 2n, N^own = 2n -I- 1 or vice versa. 

We recall once more that we are dealing with a ring, 
which is weakly connected to an electron reservoir of 
fixed chemical potential. Therefore, for non-interacting 
electrons or for an attractive electron-electron interac- 
tion between electrons carrying spin, the ring at arbi- 
trary flux will exchange pairs of particles with the reser- 
voir as the gate voltage is varied. For special values of 
the flux, when the energy levels are fourfold degenerate 
((E'/^o — integer or half-integer), this exchange will in- 
volve four electrons. Since for noninteracting electrons 
with spin an energy level is at least doubly occupied, the 
ground state (T = 0, = 0, $ = 0) in the ring must 
have an odd number of both spin- "up" and spin- "down" 
electrons. We refer to this as the * (odd,odd)-case, where 
obviously Ntot = + N°^^^^ = 4n 2. The formulas 
for the (even, even)- case (Ntot — 4n) can be derived from 
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those valid for the (odd,odd)-case by a simple shift of the 
flux by half a flux quantum (it is the parity rule for elec- 
trons with spin). The (odd,even)-case describes an odd 
number of spin-1/2 fermions in a ring (Nfot = 4n + 1 or 
4n + 3). Such a ground state can be justified only for 
a repulsive interaction. Moreover, for weakly interacting 
particles (1 — a ^ 1) it is destroyed due to the parti- 
cle exchange with electron reservoir already at very low 
temperatures Td ^ (l — a)A. 

For the (odd,odd)-case the persistent current as a func- 
tion of flux behaves very much like the current of spinless 
electrons. It oscillates with the fundamental period $o 
and amplitude Jq. At low temperatures the only differ- 
ence is the change in oscillation period when the gate 
voltage is varied. In the limit of strong repulsion the pe- 
riod of oscillations is twice of that for spinless electrons. 
The inclusion of spin changes also the high-T behavior 
of persistent current. The most striking effect is the de- 
crease by a factor of three (compared to the spinless case) 
of the crossover temperature in the limit of strong repul- 
sion. 

A ring with an odd number of spin-1/2 fermions 
(Ntot = 4n -|- 1 or 4n -|- 3) at fixed chemical potential can 
be conceived only at low temperatures. In such a situ- 
ation we find that the persistent current oscillates with 
half the period ($o/2) and half the amplitude (/o/2) com- 
pared with the case of spinless electrons. Notice that for 
non-interacting electrons the effect of period- and ampli- 
tude "halving" was predicted in Ref. |2^. We generalize 
this prediction to the more realistic case of correlated 
electrons at finite (though small) temperatures. 

The Coulomb staircase for strongly interacting spin- 
1/2 fermions {a 0) at low temperatures looks like the 
one for spinless particles. The charge in a ring is quan- 
tized in units of the elementary charge, e, and all steps 
are of equal width. The spin degrees of freedom influence 
the Coulomb blockade effects at moderate and weak in- 
teraction strengths. The width of steps corresponding 
to the number Ntot = 4n, 4ri -f 1, 4n -|- 3 of electrons in 
a ring is decreased when the interaction strength is de- 
creased and is contracted to a point in the limit a 1. 
For non-interacting spin-1/2 electrons the total charge in 
a ring at small flux values is quantized in units of 4e (as 
it should be from the qualitative considerations). 

II. THE MODEL 

In order to study Coulomb blockade effects in a ID 
ring of strongly correlated spinless electrons, we will use 
the Luttinger liquid (LL) modebllil For our purposes it is 
convenient to represent the Lagrangian of the model in 
the form: 

(1) 



Here mg is the electron (band) mass, J) = Nq/L is the 
mean density of "vacuum" electrons in the ring. [We as- 
sume that the LL-ring is in weak contact with a reservoir 
of 2D electrons. Then for repulsive interaction, Nq is 
fixed by the chemical potential /i of the reservoir which 
is implicit in the formulation Eq. (|^)]; s is the plasmon 
velocity, $ is the flux of magnetic field ($0 = hc/e is 
the flux quantum), N is the total number of particles in 
the ring, = eVg {Vg being the gate voltage). The La- 
grangian (|l|) describes the long-wavelength dynamics of 
the electron displacement field u(x,t) {(p — 2TrJ)u{x,t)). 
The total number, N, of electrons in the ring can be rep- 
resented as a sum of "vacuum electron" number, Nq, and 
the number, to, of extra electrons which appear due to 
the particle exchange between the reservoir and the ring, 
N = Nq + m. The integers m = 0,± 1, ± 2,...) can 
be connected with a winding number of topological ex- 
citations of the <y9-field in a ring geometry. It enters into 
the boundary conditions imposed on the dynamical field 
ip{x,T) in (x,r) space (r is an imaginary time) 

ip{T + (3,x) ~ (p{T,x) ^2TTn 

(^(r, X + L) ~ </7(t, x) — 27rm. ^ ' 

Here /? = 1/T is the inverse temperature, n = 0, ± 1, 
± 2, ...is the winding number for the imaginary time 
"evolution" . The twisted boundary conditions (g) are 
analogous to the ones-derived in Ref. Notice that in 
Haldane's descriptionO of the LL, the homotopy indices 
(n, to) represent topological current and charge excita- 
tions. At nonzero flux <&, the parity rule (see Ref. P) that 
connects the two sectors of the model (charge and cur- 
rent) can be explicitly incorporated into the Lagrangian 
(|l|) by adding to the external flux, <&, a fictitious (statis- 
tical) fiux = ($0/2) for an odd (even) total number 
of electrons in the ring. This statistical flux is represented 
by the second term in the parenthesis in Eq. (Q). 

For the following analysis, it is convenient to repre- 
sent the thermodynamic potential, f2, (grand canonical 
ensemble) as a path integral over the fluctuations of the 
dynamical fleld ip{x,T) 

= -Tin I 1 2?^m,«e-^-(^— I , (3) 

m.n— — OQ ) 

where Se is the Euclidean (imaginary time) action for the 
Lagrangian (|l|) and we explicitly include in the definition, 
a sum over the homotopy indices (TO,rt), which repre- 
sent the contributions of topological excitations, plasmon 
modes. 

Topological trajectories (zero modes) satisfying the 
twisted boundary conditions (^ take the form 

T X 

(Pm,nix,T) =2Tr-n + 2TT-m. (4) 

A more general solution can be obtained by adding to a 
zero mode of fixed to and n an additional part which then 
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must have periodic boundary conditions. This physicaUy 
corresponds to having an intermediate state that contains 
plasmons. Because of the periodic boundary conditions 
these extra parts only depend on temperature. Since we 
are interested in derivates of the thermodynamic poten- 
tial with respect to flux (to get the persistent current) 
and gate voltage (to get the average charge) we wilLcon- 
sider only the zero mode solutions in what follows.E3 

So, since only zero modes are important, our problem 
reduces to the evaluation of a double sum in Eq. (||). 
Since the current and charge sectors of the model are 
connected only by a single (parity) term in Eq.(|l]), this 
summation can be done exactlyliJ The result is 



47r ^ 



(5) 



93(6, q)d3{jl, q^" ) + 04(6, q)e^{% g"')] } 



Here O^^ilv.q) are the Jacobi theta- functions (see e.g. 
Ref. ^ and 



g = exp(-7r2j). (6) 



The dimcnsionless parameter a = vp/s, where vp = 
TThNo/rnQL, characterizes the correlation properties of 
the LL. 



III. PERSISTENT CURRENT 

The persistent current in a ID LL ring and its depen- 
dence on enclosed magnetic flux and gate voltage can 
be obtained from the thermodynamic potential (|^). One 
finds that 



eTJ_ 



(7) 



de 



In 



At first we consider the case of free spinless electrons 
(a=l). In this case Eq. (^ can be shown to coincide with 
a well-known expression for the persistent current of non- 
interacting electrons at fixed chemical potential that was 
derived in a standar d . ap proach using the Fermi-Dirac 
distribution function.E^O In order to show this equiv- 
alence we take advantage of the fact that for a—1 the 
sum of 0-functions in Eq. (|^) can be written as a prod- 
uct of ^-functions with the help of addition formulae for 
^-functions.tS The derivate with respect to flux can then 
be performed to give the result 



//,,e($,^5,T) 
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valid in the limit a = 1 , In Eq. (||) 6 {v, q) is the deriva- 
tive of the theta-function with respect to its argument z/, 
and we have used the reduced notation 
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(9) 



By making use of jLhe formulae for logarithmic derivatives 
of the ^-functional^ one finds 



eTA « sin(27rn-f 



Ifree - ^ 2^(-l) 
n—1 



lcos(27rn^) 



sinh(27r2n^) 



(10) 



This expresion coincides with the well-known expression 
for the persistent current of free spinless electrons at 
fixed chemical potential ^ {fj, = ep + eVg, ep is the 
Fermi energy). Notice that the crossover temperature 
T* — A/2tt'^ — Tivp/ttL in Eq. ( p^ ) is reduced by a factor 
2 in comparison with the one for an isolated ring (fixed 
number of particles). In what follows we will show that 
for interacting fermions (a ^ 1) the crossover tempera- 
ture — for a ring at fixed chemical potential — smoothly 
interpolates between the result for a stiff Wigner crystal 
(when fiuctuations of the electron number are strongly 
suppressed) and the result for a Fermi gas of spinless 
non-interacting electrons. 

Now we proceed to the case of interacting fermions. 
We will use the general result (|^) to find analytic ex- 
pressions for the low- and high temperature limits and 
rely on numerical evaluations of (|^) in the intermediate 
temperature region. At low temperatures, oscillations of 
the persistent current as a function of gate voltage (for 
$ ^ 0) are shown in Fig. |2^. The traces shown refer to 
a ring with a fixed number, A'o, of "vacuum" fermions. 
It is easy to see from the general expression Eq. (||) that 
the result for an "even" ground state can be obtained 
from those corresponding to Nq being an odd integer by 
a shift of the gate voltage, Vg ^Vg + Ajleo? . Therefore 
in what follows we will analyze only the case when the 
number of "vacuum" electrons is odd. 

The oscillations of the persistent current when the gate 
voltage is changed in the presence of a small magnetic 
flux amount to a series of transitions between a dia- and a 
paramagnetic response. At low temperatures these tran- 
sitions have the form of periodic rectangular bumps as 
shown in Fig. ||(a). The width in gate voltage of the 
bumps and their positions depend on the correlation pa- 
rameter a of the LL and the dimcnsionless flux ^/^o- 
The ratio of the width of paramagnetic response to the 
width of diamagnetic response can be simply expressed 
as 



(p) 



1 



*0 



AK 



1 



(11) 



For a stiff Wigner crystal-ring (a ^ 1) the width of 
the bumps and the dips are equal since the ground state 
of a strongly repulsive LL-ring does not depend on the 
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parity of the particle number. With less stiffness (larger 
a), the width of the current bumps (corresponding to a 
paramagnetic response) decreases and for a ^ 1 they 

reduce to narrow spikes located at vj^^ = (2k + l)A/2, 
where k is an integer. For an attractive electron-electron 
interaction (a > 1) , the oscillatory dependence on gate 
voltage vanish. The most plausible values of the correla- 
tion parameter should be a ^ 1 (we have in mind possible 
experiments in AlGaAs heterostructures). In this regime 
the shape of the current oscillations is sensitive to a and 
the mere detection of such oscillations could be a strong 
argument in favour of a non- Fermi liquid like behavior of 
quantum rings. 
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FIG. 2. Gate voltage dependence of the persistent cur- 
rent for a system of electrons without (a) and with (b) spin. 
Traces have been vertically offset for clarity. In both cases 
the electron density is the same and the gate voltage has been 
scaled by A, the level spacing at the Fermi level in the spinless 
case. If a = Oc = 1 (free electrons) one can see a doubling 
of the oscillation period in the case without spin compared to 
with spin. This is because if the densities are the same, the 
actual level spacing at the Fermi level in the spinless case is 
twice as large as in the case with spin. In contrast, for the 
case of strong interaction (a, Oc small) the ratio ac/a -^1/2 
and hence the oscillation period for electrons with spin is now 
half the period found for the case of electrons without spin. 
[An explicit expression for a in terms of the interaction po- 
tential Vb can be derived along the lines of section V. The 
result a = [1 + Vo/ttJivf]''^^^ differs from Eq. (32) for the 
case of electrons with spin only because the Fermi velocity is 
twice as large for spinless electrons of the same density]. 

In the high temperature regime T > T*{a), where 
T*(a) is a crossover temperature to be specified below, 
the flux dependence of the persistent current can be ex- 
pressed as a sum of two basic harmonics 

/((O) ^ 4^ |(_i)^0e-=(i+"=)T/A^Qg(27ra2^)x 



sin(27r^)-e-4-=^/^sin(4.^)j 



(12) 



According to Eq. ( [T^ ) , the main contribution to the cur- 
rent is given by the first term if < 3. The current 
then oscillates as a function of flux with the fundamental 
period $o — hc/e. At the same time it is parity sensitive 
and has an oscillatory dependence on the gate voltage. 
From Eq. ( p^ one finds the crossover temperature to be 
a function of the interaction parameter a: 



T*{a) 



1 A 



(13) 



This temperature attains its maximum (T* = A/tt^) for 
a stiff (a <C 1) Wigner crystal-ring, when the fluctuations 
of particle number in the ring are suppressed and the 
system can be considered as isolated. E3iij 

If Qf^ > 3, on the other hand, the second harmonic 
in Eq. ( p^ gives the main contribution to the persis- 
tent current. There is no dependence on gate volt- 
age and no parity effects in this term. The current 
is periodic in flux with a halved period ($0/2) and is 
diamagnetic. The crossover temperature saturates at 
r*(a = 73) = A/47r2, which is the lowest crossover tem- 
perature for spinless particles. The current oscillations 
at a > \/3 resemble the Aharonov-Bohm oscillations in a 
superconducting ring (period halving, absence of parity 
effect). ta This analo gy appears also when studying the 
effects of Coulomb blockade in quantum rings. 



IV. CHARGE OSCILLATIONS AND COULOMB 
BLOCKADE 

The average number of particles in a ring connected to 
a reservoir of electrons with a chemical potential /i can 
be found by making use of the general thermodynamic 
relation 



1 dn 



(14) 



Since we have an exact expression (^) for the thermo- 
dynamic potential of a quantum ring, it is straight- 
forward to find out how AA^ — N — No{fj,) depends on 
temperature, magnetic field (flux) and gate voltage. 

At low temperatures the gate voltage-dependence of 
charge CL= eAN = ef{Vg) takes the form of a Coulomb 
staircas£] with steps of even and odd heights (in units 
of e) and different widths (Fig. |^a,b). Using Eqs. (j|) and 
( |T^ ) one can derive an analytic expression for AN{Vg). 
In the limit of zero magnetic flux one finds 



AN{JI) = 271+ X 

2 ( j- j^-^K sin(27rn7I) 

^ I y-\oc sin(27rn/7) 

where k is an integer and 



(15) 

, <'p<k+-p^ 
,k = Jl^< (J, < 



2fc+l 
2 



f(l + «^),«<l (16) 
5 ,a> 1. ^ ^ 
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According to Eq. ( |15| ) the electric charge in a quantum 
ring is quantized in units of the electron charge e if the 
interaction is repulsive (a < 1; see Fig. ||b). The step 



widths corresponding to even (AfJ"') and odd {AVg'"') 
values of AN depend on the correlation parameter a and 
the dimensionless flux <i>/<i>o. The ratio of the widths is 
described by the simple equation 



^(o)^ 



(p) 



1 



AV, 



(d) 



1 



(17) 
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Recall that when deriving expression ( [T^ ) for the aver- 
age number of ring particles, the number Nq of "vac- 
uum" electrons was assumed to be odd. Therefore the 
odd(even)-step corresponds to an even(odd) total num- 
ber of electrons in a ring. 
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FIG. 3. Gate voltage dependence of the normalized 
charge for free (a) and interacting (b) electrons without spin 
at different temperatures (in units of To = A/2n, where A is 
the energy level spacing at the Fermi level). Traces have been 
vertically offset for clarity. Note that the step width in (a) 
is the sum of the (unequal) step widths of two adjacent steps 
in (b). Comparing a and b curves one can see that interac- 
tions makes the Coulomb staircase more robust with respect 
to temperature. 

From Eq. (|l^) one can see that in a stiff Wigner crystal- 
ring (a <C 1) all steps are of equal width — as it should 
be - since in this limit the ground state energy does not 
depend on the parity of the particle number. As a — > 1, 
the steps corresponding to an even total number of parti- 
cles are reduced to a point and for the LL with attractive 
interactions (a > 1) the charge in the ring is quantized 
in units of 2e. The above picture of charge quantiza- 
tion (which is obvious from the physical point of view) 
explains the oscillations between dia- and paramagnetic 
currents analyzed in the previous section. Notice that the 
ratio ( p7| ) of the step widths is sensitive to a in the region 
a < 1. This permits us to suggest the simple experiment 



of measuring the charge quantization in a quantum ring 
as a sensitive method of studying correlations in a Lut- 
tinger liquid. 

At high temperatures, T ^ T*{a), the quantization of 
charge in a quantum ring is destroyed and the Coulomb 
staircase is smeared out into a weakly modulated linear 
dependence of charge on gate voltage 



AN{V, 



T) = ^ jeK, - 47rre-4-'"'^/^ sin(47ra2^)-f 

47rT(-l)^''cos(2^^)e-^'(i+"')^/^sin(2^a2^)| . (18) 

It follows immediately that when < 1/3, the peri- 
odic modulation mainly comes from the second term of 
Eq. (^8|). It depends neither on magnetic flux nor on the 
parity of "vacuum" electrons. The crossover temperature 
for charge oscillations is equal to 



T* 



< 1/3. 



(19) 



In the limit of strong repulsive interactions (a —^ 0), 
T* — > oo, which implies that a sharp charge quantiza- 
tion is restored at any flnite temperature. 

Finally, when > 1/3 the charge oscillations are de- 
termined by the third term in Eq. ( [T^ ) rather than the 
second. The oscillation period with respect to Vg is then 
doubled and the crossover temperature decreases with 
the increase of a 



^2(1 + a2)' 



> 1/3. 



(20) 



V. INFLUENCE OF LONG-RANGE COULOMB 
INTERACTIONS ON THE THERMODYNAMICS 
OF QUANTUM RINGS 

So far it has been assumed that the interactions be- 
tween electrons in the quantum ring is short ranged. Now 
we consider long-range Coulomb forces and demonstrate 
that this interaction can be included in our scheme by 
replacing the "bare" coupling constant a by an effective 
size-dependent "interaction constant" ae//(L), 



(21) 



Here A is an ultraviolet cutoff (A <C L) which in our case 
has the physical meaning of the width of the quantum 
wire. 

Long-range Coulomb interactions can be incorporated 
in our model by adding to the local Lagrangian (^ a 
nonlocal Coulomb term 



c '(x) f d ^ 



(22) 
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This term (see e.g. Refs. ^ , p^ describes the electrostatic 
interaction of charge densities p{x) = ip {x)/2Tr localized 
in a long ID quantum wire of width A. 

In the presence of unscreened Coulomb forces, the 
equation of motion for the dynamical field ip{x,T) in 
imaginary time takes the form 




= 0, (23) 



where Sq — 2e^p/mo (p = mo/L). Since we are interested 
only in topological excitations, we will seek solutions to 
Eq.(|3|) of the form 

(P{x,t) ^2Tr^n + ij{x), n = 0,±l,... (24) 

where the spatial derivative of ip {p{x) = dxip{x)/2TT) 
obeys the integral equation 



pyx) = 7 / = Cyi)- 



^[x ~ y)2 + A2 



(25) 



Here 7 = [Aa^ /'K){e'^ /Tlvf)] C(7) is an integration con- 
stant, which will be determined in what follows by the 
requirement that the desired solution ■0(2;) has to satisfy 
the twisted boundary conditions (|^) . One readily verifies 
that 



Pi.l) 



l + 27ln(^) 



(26) 



is an approximate solution of the integral equation (25). 

It should be noted that when modelling the effects of 
long-range Coulomb interactions in a quantum ring, we 
substitute the ring geometry by a line when calculating 
the Coulomb energy, Eq.(|2^). Such a substitution can 
be justified only for a long enough ring L 3> A, when 
one can neglect numerical factors in the argument of the 
big logarithm in Eq. (^). In this case the precise val- 
ues of the integration limits (L — )■ aL, a ~ 1) in the 
Coulomb integrals Eqs. ( p3|) and ( p5| ) do not affect the 
final expression evaluated in the logarithmic approxima- 
tion. It is evident also from a physical point of view that 
in a perfect ring, a homogenous distribution of charges in 
the long wavelength limit can be inferred from symmetry 
considerations. 

By determining the integration constant C(7) from 
the requirement that the desired solution ipi^) satisfies 
the twisted boundary conditions (H), one concludes that 
long-range Coulomb forces do not deform the zero modes 
which in the ideal (impurity free) ring remain of the the 
form of Eq. (^). However, the appearance of a Coulomb 
energy in the Lagrangian of the model changes the Eu- 
cledian action for these trajecories. This effect can be 
taken into account by replacing the "bare" constant a by 
the effective size-dependent coupling constant aeff{L) of 
Eq. (|l|). 



Equation ( pl| ) for aeff{L) requires a comment. Let 
us define the correlation parameter, a, of the LL as the 
value of cteff at L = A (a = Q!e//(A)). Then Eq. ( pl| ) 
takes the typical form of a running interaction constant 
in the renormalization group sense (see e.g. Ref. ^3|). As 
L goes to infinity, aef / — > irrespective of the value of 
the "bare" coupling constant a. Hence, in the presence 
of long range Coulomb interactions, a LL-system of long 
enough length behaves like a stiiffae// <C 1) Wigner 
crystal (see the discussion in Ref. ^). In particular, the 
crossover temperature for current oscillations. 



2 „2 



hvi 



(27) 



ceases to depend on a when {e"^ / hv p) In {L / X) » 1 and 
coincidfis with the crossover temperature of isolated LL- 

rings.yy 

In experiments involving quantum wires and dots in 
the 2D electron gas of GaAlAs hctcrostructures, long- 
range Coulomb interactions can be partially screened by 
a large metallic electrode (gate) which controls the den- 
sity of charge carriers. If the distance, D (£> S> A), be- 
tween the quantum ring and the gate is smaller than 
the characteristic size of the ring. Coulomb forces are 
screened on distances of the order of D and one has to 
replace L by D in Eq. (pT|). In the experiments consid- 
ered, the distance D can be regarded as a controllable pa- 
rameter. Thus one has an interesting possibility to study 
oscillation effects in quantum rings of different stiffness 
aeffiD) < 1. 



VI. THERMODYNAMICS OF QUANTUM RINGS 
FOR SPIN-1/2 ELECTRONS 

In the previous section we studied the thermodynamic 
properties of ID mesoscopic rings of strongly interacting 
spinless electrons. Now we consider the infiuence of elec- 
tron spin on the persistent current and charge quantiza- 
tion in a quantum ring. It is well known that in a perfect 
(impurity free) Luttinger liquid charge and spin degrees 
of freedom are separated and their dynamics can be de- 
scribed by independent quadratic Lagrangians (see e.g. 
Ref. |2^). Though the local dynamics of charge and spin 
excitations is independent, globally the two sectors of the 
LL model are connected by the "parity rules" which re- 
flect the simple fact that the total spin and charge of 
the system are determined by the numbers of spin- "up" 
and spin- "down" electrons. Therefore the spin degrees of 
freedom affect persistent current- and charge oscillations 
through the parity term in the Lagrangian. 

It is easy to generalize the model Lagrangian (0) to de- 
scribe a Luttinger liquid of spin-1/2 electrons. Let N^'^j 
be the number of "up"- and "down"- spin electrons in 
the ground state (T=0, Vg=0). We will assume that the 
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electron system under consideration is not spin polar- 
ized on a macroscopic scale. It means that the densi- 
ties of electrons with opposite spin projection are equal 
Pi ~ pj^ = p/2, where p = (iv|°'' -I- iv|°^)/L is the mean 
electron density fixed by the chemical potential /i of the 
electron reservoir. If we neglect the difference in inter- 
action of electrons with different spin projections, the 
electron system in question can be regarded asi-^ two- 
component Luttingcr liquid of spinless electrons,tj 

mo 



$ N, 



h 



(0) 



nij - 1 



27r 



(28) 



Here vf — as above — refers to the Fermi velocity of 
spinless electrons of the same density. The actual Fermi 
velocity here is therefore w_f/2 = 7r?ip/2mo; U(x-y) is the 
kernel of the electrostatic electron-electron interaction. 
In what follows we will assume it to be localcj 



U{x -y) = VoS{x - y) 



(29) 



The topological excitations which determine the current 
and charge oscillations in the 2-componcnt LL arc of the 
form 



2'KT 2'KX 
ifij (X, t) = Hj — ' 



(30) 



where (nj,mj) are winding numbers, which indepen- 
dently run over integer values. The Euclidean action for 
the trajectories (30) takes the form 



{^2 
_ ... ' 
+ I/3A (1 + Wo/TThvp) + Pligrfij] 

-m^mi, (31) 



L 



where 



A=— Q. = ^ + 2 ' = (32) 

When calculating the thermo- 

dynamic potential il.{^,Vg,T) of a ID ring of strongly 
correlated electrons 



n-Qp = -T\n< ^ cxp [-5ij(ri|^l,m|^J/n] 



m-j- ^1 ,711- .| — — oo 



(33) 



it is convenient to replace the sum over (m|_|, ?if,j.) by a 
sum over linear combinations of winding numbers which 
describe the charge (c) and spin (s) channels: 



(34) 



The summation in Eq. (^3|) can be performed ecactly. We 
omit the straightforward but cumbersome intermediate 
calculations and simply give the final expression for the 
oscillating part of Q. The form of this expression depends 
significantly on the parity of the numbers (iv|"' , iYj"' ) 
of spin-" up" and "down" vacuum electrons in the ring. 
Because of this we separately analyze the cases when 
the above numbers are of equal parity - (odd,odd) and 
(even, even) or they have opposite parities (even, odd) or 
vice versa. To compare the results for electrons with and 
without spin we express all quantites in terms of Iq de- 
termined for the spinless case: 



/n = 



(35) 



VII. QUANTUM RING WITH AN ODD 
NUMBER OF ELECTRON PAIRS 

For the case of an odd number of electron pairs (spin- 
1,1), the oscillating part of the thermodynamic potential 
il takes the form 



-4:ai 



A 



(36) 



Tin- 



[ejif, q') [e^iJl, g"')03(O, q) + e.iJi, g"')04(O, q) 

+Ol{f, q^) [e^i-p, g"')04(O, q) + OiiTI, <7"')03(O, q) 

+4(?4(2/,g^)f?4(27I,9""')^?I(0,g^)}, 

where / = '^/^o is the dimensionless flux, /J = 
2a^eVo/A; the parameters q and A are defined by 
Eq. (|^. The correlation parameter Uc = vp/vc in the 
charge channel — where Vc is the velocity of charged ex- 
citations — is 



x-1 



(37) 



N^^^)/moL). The parameter ac 



where vp 

describes the correlation properties of a LL of spin- 1/2 
electrons. In the absence of any Pauli interaction of the 
electron spin with the magnetic field, the stiffness in the 
spin system coincides with that of a Fermi gas of nonin- 
teracting particles, ag = vp/vg — 1. These simple equa- 
tions for the correlation parameters of the spin- 1/2 LL 
were first derived in Ref. Ea,Eq. 



To begin with we analyze Eq. (36) in the limit of non- 
interacting electrons {Vq 0, ac 1). In this case 
Eq. (51) reduces to 
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n°!^{ac= 1) = -rin{2exp(27l) x (38) 

By making use of the addition formulae for the Jacobi 
theta-functions (see e.g. Ref. |l^) we get the follow- 
ing expression for the persistent current of free spin- 1/2 
fermions at finite temperature (Cf. Eqs. (M) and (M)) 



eT 1 



X 



h n 

2\a (— „2^ 



(39) 



- In {03(/, g')f?3(M, q^) + ^4(/, q^)94{Ti. ?')} 

= ^.„ sin(2^n^)cos(47rnS^) 
n ^ sinh(4^2,,T) 

This result coincides with the one derived by Kulik0 m 
a standard approach using the Fermi-Dirac distribution 
function /fd(p«), 



/($,/.,r) = -c5]^/FZ5(p„) 



9$' 



(40) 



er^sin(2^fc|-)cos(47rfcf) 



/c=l 



sinh(fc 



A I 



where e is the one-particle spectrum. Since the persistent 
current is a property of the Fermi surface, it does not 
depend on the precise form of the spectrum far from the 
Fermi momentum. When expressed in terms of the Fermi 
velocity, itj-is the same for quadratic and linear energy 
dispersionO 

The standard expression (^0|) for the persistent current 
in the case of noninteracting electrons with spin is trans- 
formed to our result (p9|) if we put the chemical potential 



equal io [i — -I- eVg {sp^ is the Fermi energy of the 
ground state with an odd number of spin- "up" , "down" 
electron pairs). From a physical point of view, it is ev- 
ident that for noninteracting spin-1/2 electrons, the en- 
ergy levels are doubly degenerate. Therefore the ground 
state (T — 0, Vg=0, $=0) of free electrons at fixed chem- 
ical potential has an odd number of electron pairs. We 
will prove this fact explicitly below when studying the 
effects of the Coulomb blockade. 

Notice that at T = and Vg=0, Eq. (||) describes the 
persistent current for a system with a fixed number of 
particles. In the case considered (odd number of pairs), 
our formula exactly coincides with the one obtained in 
Ref. 

We conclude that the persistent current in a ring with 
an odd number of spin- "up" - "down" pairs oscillates with 
the fundamental period $o = hc/e. It is diamagnetic and 
at low temperatures has the same amplitude Iq = evp/L 
as the current of spinless fermions (at the same particle 
density p). The influence of spin on the thermodynamic 
properties of quantum rings of free electrons is straight- 
forward — the crossover temperature and oscillation pe- 
riod on gate voltage are diminished by a factor 2 (in 



comparison with the spinless case, Eq. (|10D) which is a 
trivial consequence of the "halving" of the Fermi velocity 
for electrons with spin and same density. 

What new effects do the electron-electron correlations 
lead to? At low temperatures (T 0) the windings in 
the spatial sector are frozen (due to the suppression of 
particle fluctuations) and we get a simple expression for 
the persistent current of a LL-ring at T = of the form 



2. ,.„sin(2^n|;) 



j(r = 0) = -/„^(-l) 



n ^ " A ^ 



(41) 



Although the oscillations have the same sawtooth-like 
shape as for noninteracting particles, the period of the 
current oscillations when the gate voltage is varied is 
drastically different in the case of strongly correlated elec- 
trons. According to Eqs. ( p7| ) and ( |4l] ) one gets for strong 
repulsion current oscillations with a period which is twice 
the oscillation period for a LL-ring of spinless fermions 
(see Fig. |b). 

Electron-electron correlations significantly affect the 
Aharonov-Bohm oscillations (even for an impurity free 
ring) in the high temperature regime. In particular, the 
crossover temperature becomes a function of the corre- 
lation parameter a. Using a small-g expansion of the 
6'-functions, one can from Eq. ( p6| ) derive the high-T ex- 
pansion for the persistent current of correlated electrons. 
As in the case of spinless fermions (see Eq. (|l2|)) it is 
sufficient to keep only two harmonics in the expansion, 



Sin ("2^^)+ e-*^'^/^ Sin (^n- 



(42) 



Note that the current oscillations are determined mainly 
by the first term in ( ^ ) if < 5. The current is dia- 
magnetic and persists with fundamental period $0. The 
crossover temperature is 



T*{al) 



, < 5 
,al>5 



(43) 



The transition to the regime where the oscillationpe- 
riod is reduced to $0/2 occurs only when the interaction 
is strong and attractive, i.e. when > 5. However, 
it should be noted that unlike in the spinless case, the 
two-component model ( |2^ ) can not be directly applied 
for describing attractive interactions. It is known (see 
e.g. Ref. |2^) that in the last case backscattering pro- 
cesses become relevant and result in a gap for the spin 
excitation spectrum. The appearance of the gap rules 
out Luttinger liquid- like behavior of spin-1/2 attractive 
electrons. Therefore, in what follows we will analyze only 
the regime of repulsive interactions {ac < 1). 
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It is worth to note here that the inclusion of spin de- 
grees of freedom changes the maximum crossover tem- 
perature, which is attained in the hmit of strong repul- 
sion Us ^ I (stiff Wigner crystal) . As one can see from 
Eq. (^), for spin- 1/2 electrons this temperature is three 
times smaller than for the Wigner crystal-ring of spinless 
fermions, Eq. (|l3|). 

Now we proceed to the analysis of charge quantiza- 
tion in a LL-ring of electrons with spin. For the case in 
question (an odd number of electron pairs: A^|°^ = ivj"-* 

=2n-|-l) the analytic formulae for N{Vg) are too cumber- 
some for practical use and we will consider only numerical 
results. Some results for electrons with spin are shown 
in Fig. |. 




a;(eV/A) 



FIG. 4. Splitting of the Coulomb staircase at low tem- 
peratures and no flux due to repulsive interactions between 
spin-1/2 electrons. Traces have been vertically oflFset for clar- 
ity. In the absence of a magnetic flux the energy levels of 
the noninteracting electrons (qc = 1) are fourfold degenerate 
(with a factor of two from spin). Hence four electrons at a 
time are entering the ring each time the gate voltage Vg is such 
that a new degenerate ring level is pulled below the chemical 
potential of the reservoir. As the strength of the repulsive 
electron-electron interactions is increased (qc < 1) the non- 
interacting energy levels become increasingly irrelevant and 
electron correlations dominate. As shown this implies that 
electrons are entering the ring one by one, and in the limit 
Qc <^ 1 the steps in the Coulomb staircase become equidis- 
tant. 

In the strongly repulsive limit {ac ^ 1) and at low 
temperatures the spin degree of freedom is insignificant 
and we have the same Coulomb staircase as for spinless 
electrons. The charge quantization for the two cases con- 
sidered (electrons with and without spin) differs drasti- 
cally in the opposite limit — weakly interacting particles. 
For spin-1/2 fermions, the steps pertaining to the num- 
bers N = An, An + 1, An + 3 electrons on the ring are 
decreased with the decrease of interaction strength and 
disappear (the width of the steps is contracted to a point) 



at Qfc —> 1. 

The above picture demonstrates the physical fact that 
at fixed chemical chemical potential (grand canonical en- 
semble) the stable ground state of the Fermi gas of spin- 
1 /2 fermions always contains an odd number of spin- "up" 
and spin- "down" electrons N^^^ — A^^j"-* = 2n+l. It leads 
to a diamagnetic response of a ring when applying an ex- 
ternal magnetic field. 

Since the width of the other steps goes to zero as etc ^ 
1, one can estimate the destruction temperature of the 
corresponding "vacua" for ^ 1, Td ~ (1 — a^) A/2 
(A/2 is the level spacing for the Fermi gas of electrons 
with spin). 



VIII. QUANTUM RING WITH AN EVEN 
NUMBER OF ELECTRON PAIRS 

The thermodynamic potential, for an even num- 

ber of pairs of ring electron with opposite spin projec- 



tions N^°^ = N^"' =2n can be derived from Eq. (|36|) by 
a shift of flux, / i—rr 1/2. This is the parity rule 
for electrons with spinB3 In general, the thermodynamic 
potentials (or more correctly their oscillating parts) for 
the two cases in question — (odd, odd) and (even, even) 
— obey the following symmetry relations: 

r!;-(/,7l) = r!--(/ + 1/2,m) = (/,/!+ 1/2). (44) 

By shifting arguments in the expression ( |38| ) valid for 
an odd number of noninteracting electrons we find the 
persistent current for noninteracting electrons in a ring 
with an even number of pairs to be 

eT^ sin(27rn£)cos(4^n^) 



r(0) 



n=l 



This expression transforms 



smh{4:n'^n-^) 

(45) 

into the standard result 



Eq. ( ^ ) for ^ — ep + eVg, where ep is the Fermi en- 



ergy of the ground state with an eve n nu mber of electron 
pairs. As was discussed in section VII, such a ground 
state (A^|°^ = N^^'' =2n) for noninteracting particles is 
destroyed due to particle exchange with the reservoir. 
Thus Eq.(^) makes sense only in the T=0, Vg=0 hmit: 



free 



sin(27m^) 



(46) 



when it describes the persistent current of a ring with 
fixed {N = An) number of particles. Eq.(^) coincides 
with the one obtained in Ref. The current is para- 
magnetic, $o"P6riodic and has the same amplitude, /q 
as the persistent current of spinless fermions. The tem- 
perature behavior of the persistent current for a LL-ring 
in contact with a reservoir is described be Eqs. ( ^2|) and 
(^), independent of the parity of the number of "vac- 
uum" electrons {N\°\Nf^). 
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IX. QUANTUM RING WITH AN ODD NUMBER 
OF PARTICLES 

For a ring with an odd number of spin- 1/2 fermions 
(^jsfio) = 477, -)- 1 or 4n + 3), the oscillating part of the 
thermodynamic potential fiT'i: , takes the form 



exp 



4a 



(47) 



.3(2-,g4)04(4a^^,g4"^)- 



.4(2^,,^)03(4a?^,,^"^) 



A 

tig 

A 



The principal difference between Eqs. (^) and ( ^ ) is 
concerned with the periodicity in flux and gate voltage. 
Using the properties of the Jacobi theta-functions, one 
can immediately see from Eq. (^) that in a ring with an 
odd number of particles, the persistent current is periodic 
in magnetic flux with period <&o/2, half the fundamental 
period. 

Period "halving" for the persistent current of an odd 
number of spin- 1/2 noninteracting fermions was pre- 
dicted (for T = 0) in Rcf. Fr om the general properties 
of Eq. (47) we can conclude that this effect survives in the 
presence of interactions (see also Rcf.^O|,^l]) and even for 
finite temperatures. However, Eq. (^^ pertains to the 
case of an odd number of "ground state" electrons in the 
ring. Therefore one has to be sure that vacuum states 
corresponding to an odd number of particles {N = 4n + 1 
or 4n -|- 3) are stable against particle number fluctuations 
since we study the situation with a fixed chemical po- 
tential. We have already seen that it is not the case for 
noninteracting spin-1/2 fermions. In the last case, the 
ground state corresponds to A^'^*'^ = 472 + 2 number of 
electrons. Thus the Uc ^ 1 limit of Eq. (^ ) can be jus- 
tified only for strictly zero temperature and Vg=0 when 
it describes properties of a ring with fixed (odd) num- 
ber of noninteracting electrons. In this limit we get from 
Eq. ( ^ ) a simple expression for the persistent current. 



= 0,K, = 0) = 



1 ~ sin(47rn|^) 

-Jo > , 

TT ^ — ' n 

n=l 



(48) 



which coincides exactly with the one in Ref. The cur- 
rent is paramagnetic and periodic with half the period, 
<I>o/2, and half the amplitude, i/p, in comparison with 
the case of spinless electrons. 

For repulsive interactions, the "odd" ground state is 
stable only at temperatures T ^ ~ (1 — Q!^)A. There- 
fore it is reasonable to study current oscillations for the 
(even-odd)-case only in the low temperature regime. The 
high temperature behavior of the current is always de- 
scribed by Eq. (^). 

At low temperatures the oscillations of current with 
varying gate voltage at <I> ^ can be interpreted as os- 
cillations between a dia- and a paramagnetic response. 



Qualitatively these oscillations are of the same form as 
for the spinless particles (see Fig. ||). For strong repul- 
sion (ckc ^ 1 , stiff Wigner crystal) the bumps of param- 
agnetic response and dips of diamagnetic response have 
the same width in gate voltage. This property is closely 
related to the effects of charge quantazation. For a stiff 
Wigner crystal, the "even" and "odd" steps of the func- 
tion eAN{Vg, T ^ 0) = -{dn+-/dVg) have equal widths 
since in this strongly repulsive limit considered the en- 
ergy does not depend on the parity of the electron num- 
ber. With the increase of ac the steps corresponding to 
an odd total number of electrons on a ring get narrow 
and paramagnetic bursts of current reduce to spikes. In 
the region Qc ^ 1 (even,odd), the ground state becomes 
absolutely unstable and to study the persistent current of 
noninteracting electrons at fixed chemical potential, one 
has to start from a different expression for Q, namely 
Eq. (^), "built" on the stable ground state. 
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